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The  Seismo-Acoustic  Disturbance 
Produced  by  a  Titan  lll-D 
With  Application  to  the  Space  Transportation 
System  Launch  Environment  at  Vandenberg  AFB 


I.  INTRODUCTION 

1.1  Purport 

This  study  was  conducted  In  an  effort  to  forecast  the  seismic  and  acoustic 
fields  for  Space  Transportation  System  (STS)  launches  at  Vandenberg  Air  Force 
Base  (VAFB)  In  a  form  that  permits  vibration  estimates  within  critical  STS 
ground  facilities  (abbreviations  used  in  this  report  are  contained  In  the  List  of 
Abbreviations  following  the  text). 

1.2  Approach 

With  SAMSO  support,  AFGL  undertook  the  following  Items  of  work  to  develop 
a  technology  base  for  predicting  motion  due  to  STS  operations  at  Space  Launch 
Complex  6  (SLC  6),  VAFB.  1 

1.  From  physical  considerations.  Identify  those  properties  of  the  lnunch 
pressure  field  needed  to  predict  ground  motion. 

2.  Measure  the  acoustic  load  developed  by  a  large  class  chemical  rocket  in 
a  manner  needed  for  making  ground  motion  and  facility  vibration  forecasts. 


(Received  for  publication  14  November  1980) 

1.  SAMSO/AFGL  Memorandum  of  Agreement,  dated  31  January  1979. 
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3.  Forecast  those  properties  of  the  acoustic  load  due  to  a  STS  launch 
needed  for  ground  motion  and  facility  vibration  estimates. 

4.  Determine  the  ground  response  at  SLC-G  to  a  coherent  acoustic  load 
over  a  path  of  interest. 

5.  Predict  the  ground  motion  in  the  general  area  of  the  Launch  Control 
Center  (LCC)  for  a  STS  launch  at  SLC-fi. 

1.3  Summary  of  Finding* 

Our  observations  indicate  that  the  surface  pressure  developed  during  a  Titan 
1II-D  launch  travels  largely  as  a  nondispersive,  spacially  coherent  acoustic  load 
over  the  distances  and  frequencies  of  concern.  Also,  the  "accepted"  directivity 
index,  based  on  extrapolating  pressure  levels,  differs  markedly  from  the  plume 
radiation  pattern  based  on  phase  velocity.  This  study  also  determined  that  the 
motion  of  the  ground,  excited  by  a  rocket  launch,  is  dominated  by  the  complex 
branch  of  the  so-called  "air-coupled"  wave. 

An  admittance  term  obtained  from  small  atmospheric  explosions  is  used  to 
predict  motion  at  S.LC-U  for  a  STS  launch.  The  predicted  motion  is  sensitive  to 
our  forecast  of  the  pressure  level,  phase  velocity  and  spacial  coherency  of  the 
acoustic  load  generated  by  the  STS  plume  over  the  path,  time  and  frequencies  of 
interest. 


2.  PHYSICAL  CONSIDERATIONS 
2.1  Empirical 

Ground  motion  excited  by  a  rocket  launch  travels  outward  along  one  of  two 
paths.  The  first  path  runs  directly  through  the  ground  Trom  the  stressed  zone  in 
the  immediate  vicinity  of  the  launch  pad.  The  second  and  more  important  path 
is  determined  by  coupling  the  acoustic  load  into  the  ground  along  the  travel  path. 

Commonly,  the  large  density  contrast  between  the  air  and  ground  effectively 
decouples  the  acoustic  and  seismic  regimes.  The  normal  acoustic  impedance  Z, 
defined  by  the  amplitude  ratio  of  the  surface  pressure  to  the  ground's  vertical 
particle  velocity,  typically  is  very  large.  ^ 

Z(c,  f.  7.  0)  -  pic,  f)/V(c,  f)  LARGE 


2.  Morse,  P.  ,  and  Feshbach,  H.  (1953)  Methods  of  Theoretical  Physics.  Part  II. 
McGraw-Hill,  New  York.  “  .  * 
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Conversely,  the  reciprocal  of  impedance,  admittance,  is  quite  small. 


Z(c,  f,  z  -  Of1  -  Vic.  f)  -  V/p 

For  the  usual  condition,  the  ground  behaves  nearly  as  a  rigid  body  to  acous¬ 
tic  inputs.  The  incident  pressure  p(o,  f)  is  completely  reflected  without  a  phase 
change,  and  the  surface  pressure,  p^ie,  f,  z  =  0)  is  double  the  incident  pressure. 

Ps<<\  f.  0)  =  pj  ^  P,,  -  2pt 

Measurements  of  rocket  launches  at  the  NASA  Kennedy  Space  Center  by 

•i  4 

Mickey'  and  McCarty  and  Daiins  show  the  seismic  disturbance  coming  directly 
through  the  ground  from  the  launch  zone  to  be  negligible  compared  to  the  seismics 
following  the  acoustic  arrival.  We  also  find  this  to  be  true  for  our  observations 
at  VA I  B;  ground  motion  prior  to  the  acoustic  arrival,  carried  directly  by  the 
high  speed  seismic  ground  path,  is  small.  As  at  the  Kennedy  Space  Center, 
rocket  launch  generated  seismics  are  due  to  acoustically  excited  ground  reson¬ 
ances  whose  characteristics  are  deto  mined  by  the  elastte  properties,  densities, 
and  layering  of  the  geological  section  surrounding  the  launch  area.  The  ground 
structure  both  at  SLC-4  and  SLC-(i  is  also  of  a  kind  that  supports  "air-coupled" 
seismic  waves. 


2.2  'iVort-tieal 

2.  2.  1  PLANK  WAVL  ACOUSTIC  EXCITATION  OF  A  LAV  EH 

Insights  auout  the  relation  between  surface  pressure  and  ground  motion  can 
be  developed  by  considering  the  motion  of  a  low  speed  "alluvial"  layer  sandwiched 
between  " atmospheric"  and  "rock"  half-spaces  in  the  manner  of  Figure  1.  As 
shown,  we  excited  the  model  from  above  by  a  homogeneous  acoustic  plane  wave 
of  frequency  f,  phase  velocity  c,  and  pressure  amplitude  p. 

C, 

Modifying  the  formulation  of  llaskell  and  the  computational  algorithm  of 
Abo-Zena'1  for  a  layered  half-space  (Appendix  A),  we  compute  the  magnitude  of 


3.  Mickey,  W.  (1502)  Seismic  Investigation  of  Missile  Launchings.  U.S.  Dept. 

of  Commerce,  Coast  and  Geodetic  Survey  Publication. 

4.  McCarty,  V.,  and  Dallns,  I.  (1874)  Frequency  shift  in  air-coupled  surface 

waves  originated  by  rocket  launches,  J,  Geophys.  Res.  7 0 ; 7 0 2 7 -7039. 

5.  Haskell,  N.  (19(12)  Crustal  reflection  of  P  and  SV  waves,  Bull,  Seism.  Soc. 

Am.  £7:4751-4797. 

9.  Abo-Zena,  A.  (1979)  Dispersion  function  computations  for  unlimited  frequency 
values,  Gecphys.  J.R.  Astr.  Soc.  58:91-105. 
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Figure  1.  Alluvial  Layer  Model 


the  admittance  at  the  "air-alluvial"  interface  for  an  acoustic  input  near  grazing 
incidence  c  «  c  as  seen  in  Figure  2.  We  choose  to  compute  admittance  here, 
for  it  is  in  harmony  with  our  concept  of  the  response  of  a  linear  elastic  system 
to  an  "input  pressure"  that  excites  an  "output  particle  velocity,  '  In  Figure  2,  the 
magnitude  of  admittance,  the  amplitude  ratio  of  the  vertical  particle  velocity  to 
surface  pressure,  is  given  as  a  function  of  frequency  for  the  model.  As  depicted, 
the  magnitude  is  quite  small  for  all  but  a  subset  of  frequencies,  the  so-called 

*7 

"air-coupled"  frequencies.  The  value  of  admittance  and  the  location  of  the  air- 

coupled  frequency  is  sensitive  to  the  elastic  parameters,  layer  thickness  as  well 

7  8 

as  the  phase  velocity  of  the  acoustic  load.  ’ 

We  now  locate  large  admittance  values  of  our  model  as  they  relate  to  the 

phase  velocity  and  frequency  of  the  acoustic  input.  The  loci  of  high  response 

values  are  found  to  cluster  along  various  modal  lines  under  this  mapping.  For 

phase  velocities  cR  s  c  <  c  the  loci  or  high  responses  coincide  with  the  locked 

normal  modes;  that  is,  the  inhomogeneous  Rayleigh  waves  characteristic  of  a 
V 

free  surface  layer.  For  c  >  c  ,  the  alluvial  layer  behaves  like  a  leaky  wave 

7  q  °  ' 

guide.  ’ 


7.  Ewing,  W.,  et  al  (1957)  Elastic  Waves  in  a  Layered  Media,  McGraw-Hill, 

New  York.  “  — — — 

8.  Crowley,  F. ,  and  Oasing,  H.  (1969)  On  the  Application  of  Air-Coupled 

Seismic  Waves.  AFC'RL-69-03 12. 

9.  Gilbert,  F.  (1964)  Propagation  of  transient  leading  modes  in  a  stratified 

elastic  waveguide,  Reviews  of  Geophysics  3:123-153. 


asss 
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Figure  2.  Admittance  Magnitude:  Air  Alluvial  Interface  Model 


A  a  caj  be  aeen  In  Figure  3,  large  admittance  values  generally  migrate  toward 
smaller  frequency  values  for  Increasing  phase  velocity.  Looking  anead  to  the 
ground  motion  excited  by  a  rocket  launch,  we  can  anticipate  that  the  frequency  of 
air-coupling  will  also  migrate  as  the  phase  velocity  of  the  Incident  acoustic  wave 
excited  by  the  rocket  changes  with  time. 

Attributes  of  the  air-coupled  phase  have  been  discussed  at  length  in  the  seis¬ 
mic  literature.  ^  In  particular,  the  phenomena  of  surface  pressure  at 

12  13 

Rayleigh  velocities  has  been  considered  both  by  Haskell  and  Gubanov.  To 

show  the  relation,  we  now  compute  the  ratio  of  surface  pressure  to  Incident 

pressure  for  our  model  for  c  «  e  .  Here  we  find  the  computed  value  of  surface 

<> 

pressure  to  be  quite  close  to  double  the  incident  source  pressure  almost 


10.  Donn,  W. ,  et  al  (1971)  Air-coupled  seismic  waves  at  long  range  from  Apollo 

launchings,  J.  R,  Astr.  Soc.  20:101-171. 

11.  Espirosa,  A.,  and  Mickey,  W.  (1908)  Observations  of  coupled  seismic  waves 

from  sonic  booms,  a  short  note,  Aeustica  20:88-91. 

12.  Haskell,  N.  (1951)  A  note  on  air-coupied  surface  waves.  Bull.  Seism.  Soc. 

Am.  41:295-300.  - 

13.  Gubanov,  A.  (1945)  Rayleigh  waves  on  a  boundary  between  a  solid  and  a 

liquid,  J.  Exptl.  Theoret.  Phys.  (USSR)  15:497. 
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Figure  3.  Layer  Response 


everywhere  Ir.  the  manner  of  a  rigid  earth,  except  at  the  alr-eoupled  frequencies 
(see  Figure  4).  At  the  ulr-coupled  frequencies,  the  reflected  acoustic  wave  Is  * 
out  of  phase  with  the  Incident  term;  the  surface  pressure  approaches  aero.  In 
our  model  we  find  that  the  relation  between  surface  pressure  and  free  field  pres¬ 
sure  Is  a  nenatllve  function  of  the  model's  geologic  parameters.  High  resolution 
surface  acoustic  spectra  also  should  be  site  sensitive. 

2.  2.  2  CONCENTRATED  SOURCE 

The  phenomenon  of  the  air-coupled  phase  excited  by  a  concentrated  pressure 

7 

source  has  been  considered  by  Ewing  et  al  as  an  extension  of  Rayleigh  wave 
propagation  at  a  layered  half-space.  After  Ewing  (see  Figure  5),  we  find  the 
fundamental  ai  ’-coupled  phase  trails  the  sound  wave  as  an  almost  constant  fre¬ 
quency  vibration  with  a  duration  T  given  by 

T  -  M 1  /  >.'0  -  l/l<0> 

where  c  ,  1’0,  are  the  phase  and  group  velocities  of  the  half-space  model  at  the 
air-coupled  frequency.  As  shown  in  Figure  5,  the  dotted  line  depicts  the  departure 
in  the  normal  mode  response  of  the  half-space  model  due  to  the  atmosphere.  For 
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FREQUENCY  (H>) 


Figure  4. 


Ratio  of  Surface  Pressure  to  Incident  Pressure 


Figure  5.  Air-Coupled  Rayleigh  Wave 


phase  velocities  around  the  speed  of  sound  in  air,  the  layer  response  couples  al 
the  half-space  Rayleigh  wave.  For  phase  velocities  above  c  ,  the  modal  response 
is  diminished  by  leakage.  Since  bandwidth  Inversely  depends  on  duration,  we  can 
also  anticipate  that  the  bandwidth  of  the  air-coupled  phase  will  be  Inversely  re¬ 
lated  to  the  distance  to  the  source.  In  practice  we  find  that  the  bandwidth,  or 
equivalently  the  duration  of  the  air-coupled  phase,  depends  more  on  the  uniformity 

Q 

and  extent  of  the  layering  neighboring  the  observer. 

2.3  Prediction 

From  both  model  calculations  and  field  measurements,  we  find  that  a  surflclal 
layer  having  the  density  and  elastic  constants  typical  of  unconsolidated  sediments 
substantially  alters  the  "usual"  dense  earth  relation  between  the  free  field  pres¬ 
sure  and  surface  pressure  as  well  as  the  relation  between  surface  pressure  and 
ground  motion.  The  ground  response  due  to  a  propagating  acoustic  load  can  be 
quite  intense  around  the  "air -coupled"  frequencies.  For  the  usual  condition  of  a 
low  velocity  unconsolidated  layer  overlying  a  high  speed  rock,  the  frequency  of 
the  "alr-coupled"  phase  tends  to  migrate  to  lower  frequencies  and  smaller  absolute 
values  at  higher  phase  velocities. 

Given  an  Input  acoustic  pressure,  the  output  ground  motion  ran  be  forecast 
through  admittance.  For  a  uniform  area,  the  admittance  can  be  computed  once 
the  elastic  properties,  densities,  and  layer  thicknesses  are  known.  For  an  area 
of  unknown  or  complicated  structure,  the  admittance  over  any  given  path  can  be 
measured.  We  take  the  latter  approach,  namely  premeasure  the  admittance  over 
a  path  of  Interest. 

In  order  to  forecast  seismic  motion  due  to  an  STS  luunch  at  Sl.CJ-G,  VAFli, 
through  ground  admittance,  we  must  first  establish  the  Intensity,  frequency, 
phase  velocity,  and  spaclal  coherency  of  the  acoustic  load  expected  for  a  STS 
launch.  We  attempt  this  by  extrapolating  measured  values  found  for  a  large  class 
chemical  rocket  (Titan  III -D)  to  corresponding  terms  governing  the  acoustic  load 
expected  for  a  STS  launch. 


3.  MEASUREMENTS 
3.1  Syitem  Ret  pome 

The  location  of  the  AFGL  sensors  used  to  characterize  the  acoustic  load  gen¬ 
erated  by  a  Titan  III-D  rocket  is  given  in  Figure  G.  The  sensor  array  consists  of 
nine  Low  range  pressure  transducers  and  seven  seismometers.  Factory  calibra¬ 
tions  for  the  pressure  sensors  were  certified  at  AFGL;  the  seismic  elements 
were  calibrated  in  place  at  the  launch  complex.  The  array  was  located  along  a 


14 


V  KK  IICAI 


SLC-<*E 


A  SEISMIC  SENSOR 
A- SEISMIC  SENSORS  iiomzuniai. 
I 


•  PRESSURE  SENSOR 


VIIIHO.  ACOHS  lie  AMHAV 


Figure  G.  Vlbro-Acoustlc  Array 


strip  northwest  of  the  SLC-4E  launch  pad.  The  array's  dimensions  and  location 
were  set  by  the  distance  to  and  size  of  critical  STS  support  facilities  at  SLC-6. 
The  pressure  sensors  were  located  on  the  surface;  the  seismometers  wore  burled 
In  shallow  holes. 

The  output  of  each  of  the  IG  sensors  was  amplified  and  transmitted  to  a  data 
acquisition  system  located  In  the  Launch  Operations  Building  (LOB).  Here  the 
observations  were  digitized  at  a  100  samples  per  second  rate  and  stored  on 
magnetic  tape  fo'  .uter  retrieval  and  analysis.  The  period  of  measurement 
started  1  min  before  Ignition  and  continued  to  20  min. 

The  design  response  of  a  pressure  channel  is  shown  In  Figure  7;  the  corre¬ 
sponding  response  for  seismic  measurements  Is  given  in  Figure  8. 

3.2  Protection  Ratio# 

Since  the  data  were  digitized  at  a  100  sample  rate,  our  measurement  band  Is 
limited  to  something  less  than  50  Hz.  An  estimate  of  the  useful  measurement 
band  can  be  calculated  from  the  ambient  noise  Tigure  and  the  amplitude  protection 
ratio. 

The  first  fold  protection  ratio  for  the  pressure  data,  defined  by  the  ratio  of 
the  amplitude  response  at  the  measurement  frequency  to  the  corresponding 
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Figure  7.  Pressure  Channel  Amplitude  Response 


Figure  8.  Seismic  Channel  Amplitude  Resoonse 
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response  at  the  first  fold  alias  of  that  frequency,  is  given  In  Figure  9.  The 
design  protection  ratios  for  the  seismic  and  pressure  channels  are  quite  similar. 
The  reason  for  this  similarity  is  elementary:  the  signal  and  noise  characteristics 
for  both  pressure  and  seismic  measurements  above  nur  midband  frequency  are 
similar.  Above  25  lix  both  the  pressure  and  seismic  inputs  rolloff  at  6  db  per 
octave.  For  a  midband  measurement,  we  anticipate  that  the  "effective"  protec¬ 
tion  ratio,  determined  by  the  -atio  of  the  actual  value  of  the  measurement  Inband 
to  the  corresponding  aliared  term,  is  about  three  times  better  than  that  shown  in 
Figure  9. 
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Figure  ».  First  Fold  Protection  Hatlo:  Pressure  Channels 


3.5  Mwuurement  IN  out* 

Our  estimate  of  measurement  noise  is  based  on  seismic  and  pressure  obser¬ 
vations  taken  1  min  before  ignition.  We  treat  these  data  to  be  from  a  stationary 
random  process  that  is  uncorrelated  and  additive  to  the  launch  signal.  Under 
these  assumptions  the  noise  and  signal  spectra  will  simply  add. 


14.  Blackman,  H.  .  and  Tukey,  W.  (1959)  The  Measurement  of  Power  Spectra, 
Dover  Publications,  New  York.  “  * 
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The  ambient  noise  spectral  estimates  given  here  are  based  on  the  Discrete 
Fourier  Transform  (DFT's)  averaged  over  twenty  2.  5(i-sec  measurement  Inter¬ 
vals.  The  variance  of  the  averaged  unweighted  periodogram  estimates  is  about 

15 

5  percent  of  the  square  of  the  spectra  of  the  idealized  process.  Hanued  meas¬ 
urement  noise  spectral  estimates,  referred  to  the  input,  are  given  in  Figure  10 
and  Figure  11.  The  variance  of  these  estimates  is  expected  to  be  about  fh  per¬ 
cent  of  that  of  the  unweighted  form. 

The  ambient  noise  value  given  here  is  a  composite  of  ambient  pressure  or 

ground  motion  as  well  as  amplifier,  cable  and  digitization  noise.  In  our  passband 

_2 

<1-30  Ilz),  the  HMS  total  noise  figure  for  the  seismic  channels  is  1.25  X  10  mm/ 
secj  the  corresponding  term  for  the  pressure  measurements  Is  li.  5  X  10"^  pal 
(1.0  X  10  kg/cm  ).  It  should  also  be  recognized  that  the  low  frequency  limb  of 
the  pressure  spectra  is  very  sensitive  to  wind  speed  at  the  time  of  measurement. 
After  Kimball,  1()  we  can  anticipate  a  50-fold  Increase  in  the  surface  oreasure 
spectra  for  a  2-fold  increase  in  wind  speed.  A  nom><'cuStlc  pressure  term,  due 
to  launch  induced  turbulence,  can  be  anticipated  at  low  frequencies. 


Figure  10.  Noise  Spectral  Estimate:  Pressure 


Tile  references  cited  above  will  not  be  listed  here.  See  Heferences,  page  55, 
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Figure  11.  Noise  Spectral  Estimate:  Seismic 


4.  OBSERVATIONS 
4.1  lntruductioa 

The  surface  pressure,  as  well  as  the  vertical,  radial,  and  transverse  seis¬ 
mic  motion  caused  by  the  Titan  III -O  launch  of  l(i  March  1970,  is  shown  In 
Figure  12  for  an  observation  point  890  ft  (271  m)  northwest  of  SLC-4E.  The  data 
shown  here  are  corrected  for  our  system  response  in  that  the  effect  of  the  system 
response  on  the  observations  is  suppressed  by  linear  deconvolution.  The  inver¬ 
sion  operators  used  here  are  extended  until  they  recover  the  expected  input  with  a 
mean  square  error  of  less  than  0. 1  percent  (Appendix  B). 

As  can  be  seen,  the  duration  of  the  significant  disturbance  is  less  than  half  a 
minute.  The  launch  signals  start  abruptly  with  an  "ignition  pulse"  that 's  fol¬ 
lowed  by  a  butldup  in  level  that  reaches  a  maximum  intensity  in  about  9  sec. 
Following  this  maximum,  the  plume  generated  surface  pressure  level  trails  off 
in  exponential  fashion.  As  at  Cape  Kennedy,  we  find  that  conspicuous  ground 
motion  follows  the  acoustic  arrival.  The  envelope  of  the  seismic  motion  roughly 
mirrors  that  of  the  pressure.  Also,  the  motion  intensity  level  is  pretty  much  the 
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Figure  12.  Observed  Component  Motions:  Titan  111  -  U  Source 


same  on  each  ground  component.  This  unexpected  result  is  in  harmony  with  like 
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ground  measurements  at  Kennedy. 

4.2  Ignition  Pulw 

The  observed  ignition  pulse  can  be  represented  as  u  propagating  wavelet  of 
something  less  than  2  sec  duration  (Figure  13)  with  an  attenuation  factor  (A./A^)'" 
independent  of  wave  number,  and  a  propagation  delay,  (A.j  -  A())/r  independent  of 
frequency: 

p'(A.,  t)  =  (A./Ao>'n  p(Ao,  tQ  -  (A.,  -  Aq)/c)  .  (1) 


17.  JHallard,  ft.  (1970)  Dynamic  Behavior  of  Launch  Facility  Foundation  and 

Surrounding  Areas,  Miscellaneous  Paper  S-frO-Hfl,  U,  S.  Army  Water w ays 
Experimental  Station. 
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Figure  13.  Ignition  Pulse  Wavelet 
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A  best  estimate  of  uttenuutlon  and  delay  I(Aj/A())  a(Aj  -  A  ^/e]  is  defined  to 
be  the  (n,  c)  pair  that  minimizes  the  sum  square  difference  between  the  measured 
and  predicted  pressure  values,  namely 

n 

1/n-l  £  (p(Aj,  t . )  -  p't^.t,))2  <ef>mln  .  (2) 

l-l 

In  turn,  a  figure  of  merit  M  is  formed  to  test  the  representation: 

Mj*  (*-<«J>/sJl  • 

where  S2  is  the  mean  square  of  the  measured  pressure  wavelet  found  at  A  tn 
o  1  o 

excess  of  the  uddltive  ambient  noise  existing  just  prior  to  launch. 
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4.  2.  1  VELOCITY  ESTIMATE 

Figure  14  Is  a  plot  of  the  best  delay  estimates,  as  they  relate  to  separation 

(A.  -  A()).  The  phase  delay  times  computed  here  are  given  to  a  precision  that  is 

one  order  of  magnitude  higher  than  0.01  sec  of  the  field  measurements.  This  is 

1 8 

done  by  taking  advantage  of  the  shift  properties  of  the  Fast  Fourier  Transform, 
namely 

p(Aj,t  -  T.)  =  FFT,1(FFT(p(Aj,  t)  .  el  ju,T)l  . 

The  velocity  of  propagation  fox-  the  ignition  wavelet  based  on  oui'  r  estimates  is 
1102  ft/sec  (335.9  m/soc),  the  speed  of  sound  in  air  for  the  wind  and  temperature 
conditions  prevailing  at  launch  time, 

4.2.2  ATTENUATION  ESTIMATE 

The  log  of  best  attenuation  factors  is  given  in  Figure  15  as  they  relate  to  the 
source  distance.  The  attenuation  exponent  found  here  lies  midway  between  the 
attenuation  of  a  point  and  line  source  in  a  homogeneous  whole  space.  ^  Clearly, 
the  attenuation  is  significantly  less  than  the  (1/A)  term  often  used  to  adjust 
pressure  estimates.  *‘1*  ^ 

Figure  1  L>  summarizes  the  residual  error  for  our  representation  as  it  relates 
to  separation.  The  mismatch  between  the  observed  and  predicted  ignition  wave¬ 
let  pressures  leaves  a  residual  error  of  slightly  more  than  1  percent,  our  rep¬ 
resentation  of  the  ignition  pressure  wavelet,  as  given  In  Eq.  (1),  is  found  to  be 
quite  insensitive  to  separation.  Over  the  range  of  our  measurements  the  ignition 
surface  pressure  wavelet  propugutes  coherently  and  nondlspcrsively  (without  a 
change  in  form)  at  the  speed  of  sound  in  air  with  an  attenuation  factor  that  is 
Somewhat  less  than  that  ascrlbable  to  a  spherical  divergent  wavefront. 

4.3  Admittance 

4.3.1  MEASUREMENT 

We  now  use  the  ignition  wavelet  to  establish  the  admittance  for  a  near  surface 
pressure  source  exciting  the  ground  over  our  measurement  path  (see  Figure  17): 

Y(ca,  f)  . 

18.  Oppenheim,  A.,  and  Schafer,  R.  (1975)  Digital  Signal  Processing,  Prentice- 

Hall,  N.J.  - - - “ 

19.  Cole,  J.  et  al  (1957)  Noise  Radiation  from  Fourteen  Types  of  Rockets  in  the 

1.  000  to  130.  000  Pounds  Thrust  Range.  WADC  Technical  Report  57-3  54. 

teiwf; —  - 

20.  NASA  Report  Sp-8072  (1971)  Acoustic  Loads  Generated  by  the  Propulsion 

System .  ” 
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Figure  16.  Bandpass  RMS  Surface  Pressure 
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Figure  17.  Observed  Ground  Admittance 


The  admittance  computed  here  is  the  ratio  of  the  FFT  of  the  ground  particle 
velocity  and  the  ignition  wavelet  {Appendix  C).  As  anticipated  from  theoretical 
considerations,  the  measured  admittance  is  a  sensitive  function  of  frequency. 

Over  this  path  at  SLC-4,  the  fundamental  air-coupled  mode  for  acoustic  loads 
propagating  near  grazing  incidence  (c  ~  o  ),  occurs  around  5  Hz.  A  higher  mode 
response  is  also  evident  around  19  Hz. 

4.3.2  PREDICTION 

We  now  predict  the  entire  ground  motion  sequence  due  to  the  Titan  III-D 
launch  using  a  linear  operator  that  satisfies  the  ignition  transient's  admittance. 
Figure  18  depicts  the  true  and  predicted  ground  motion  for  the  Titan  III-D  launch 
under  the  hypothesis  that  the  admittance  value  for  c  =  c^  is  adequate  to  predict 
ground  motion  generated  by  the  Titan  at  SLC-4  for  the  complete  launch  sequence. 
We  find  that  the  predicted  motion  overestimates  the  true  during  and  following  the 
time  of  peak  loading  when  the  rocket  is  at  altitude. 

Figure  19  shows  the  growth  in  the  relative  ground  motion  prediction  error 
as  it  relates  to  the  Titan's  altitude  and  launch  time.  Ground  motion  estimates 
that  are  in  harmony  with  Y(c^,  f)  produce  progressively  larger  relative  error  as 
the  Titan  III-D  gains  altitude,  particularly  above  500  ft  (152  m).  This  finding 
agrees  with  our  model  response  at  higher  phase  velocities.  The  seismic  re¬ 
sponse  found  here  is  smaller  and  shifted  in  frequency  at  the  higher  phase  velocities. 
As  in  the  model,  we  find  that  the  seismic  energy  leaks  out  of  the  surficial  layer. 

In  turn,  an  accurate  forecast  of  ground  motion  around  and  after  peak  loading 
requires  that  the  admittance  prediction  operator  Y(c,  0  reflect  the  change  in  the 
phase  velocity  of  acoustic  load  produced  by  the  rocket  at  altitude;  the  phase 
velocity  of  the  acoustic  load  generated  by  the  plume  is  an  essential  parameter  for 
making  ground  motion  forecasts.  The  same  can  also  be  said  of  elastic  structures, 
namely  that  the  phasing  of  the  acoustic  load  is  an  essential  parameter  for  making 
structural  response  estimates.21 

4.4  Hume  Emitted  Acouotici 

4.  4.  1  TIME  AVERAGING 

Figure  18  depicts  the  surface  pressure  over  the  first  30  sec  of  flight  as 
measured  at  a  station  located  some  890  ft  (271  m)  from  the  SLC-4E  launch  pad. 
Following  the  ignition  pulse,  the  pressure  level  builds  for  about  ih  sec.  By  this 
time  the  Titan  III-D  has  climbed  tc  an  altitude  of  about  1000  ft  (304.8  m),  almost 
directly  over  the  launch  point  (see  Table  D1  in  Appendix  D).  At  this  altitude  we 
assume  that  the  plume  is  close  to  steady  state  and  undeflected;  hence  the  pressure 
field  is  insensitive  to  launch  pad  geometry. 

21.  NASA  Report  GP-1059,  Revision  A  (1970)  Environment  and  Test  Specification 
Levels  Ground  Support  Equipment  for  Space  Shuttle  Launch  Complex~39:TT 
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Figure  18.  Predicted/Measured  Ground  Motions:  Titan  111 -ID 
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Figure  19.  Relative  Ground  Motion  Prediction  Error 


For  a  fixed  observer,  the  surface  pressure,  p^it),  emitted  by  the  plume  - — 

appears  as  a  transient,  existing  only  in  the  Interval  (0,  T).  In  turn,  all  statistical 

time  averages  using  p,(t),  such  as  autoconvariance,  and  consequentially  power 
s  22 

spectral  density,  are  formally  zero  ;  that  is, 


c 


PP 


(r)  =  1/T 
T  -»oo 


+T/2 

J  iJs(t)Ps(t  +  dt 

-T/2 


0  . 


A  correct  estimate  of  the  power  spectral  density  of  the  pressure  generated 
by  the  undeflected  steady  state  plume  requires  that  the  temporal  effects  due  to 
rocket  mot. on  be  suppressed.  To  do  this,  we  must  locate  and  compensate  for  the 
motion  of  the  plume  acoustic  source  in  the  observer's  coordinates. 

To  demonstrate  the  point  about  time  averages,  we  compute  an  ensemble  of 
root  mean  square  pressures  using  the  measured  surface  pressure  and  different 
averaging  times.  The  resulting  values  are  plotted  in  Figure  20  for  a  common 
start  time.  As  expected,  the  calculated  values  are  sensitive  both  to  epoch  and 
duration.  Spectral  estimates  obtained  by  time  averaging  are  affected  in  precisely 
the  same  manner.  Discordant  power  spectral  estimates  can  be  expected  when 
averaging  times  are  not  specified.  Unfortunately,  much  of  the  rocket  literature 
suffers  from  this  elementary  defect. 

4.4.2  PLUME  EMITTED  SURFACE  PRESSURE 

During  the  early  segment  of  flight,  we  represent  the  pressure  pg(t),  observed 
at  a  fixed  ground  station  due  to  the  undeflected  steady  state  plume,  to  be  the 
resultant  of  a  purely  path  dependent  function,  G(r,  tjr^,  M,  and  a  directionally 
sensitive  source  shaping  operator,  P,(t  ),  that  convolves  with  a  zero  mean, 

V  o  g 

stationary,  independent,  normal  process,  N(o,  c  ): 


po(0  =  H(t)*V^(t)*Q(r,  t/ro,  t^P^t^a  •  e(tQ)  (3) 

where  (e(to))  =  0  and  =  1  •  6 <t). 

As  represented  here,  the  plume  is  replaced  by  an  equivalent  axial-symmetric 
source*  (see  Figure  21),  with  R  the  observer-source  distance  and  $  the  observer- 
source  elevation  angle.  In  turn,  t  is  the  observer's  time  to  include  the  propagation 
delay  and  V^(t)  is  an  operator  that  modifies  the  free  field  pressure  to  account  for 
the  atmospheric-ground  structure.  H(t)  is  the  impulse  response  of  the  measuring 
system. 

-  \ 

22.  Lee,  Y.  ( I960)  Statistical  Theory  of  Communication.  John  Wilev  and  Sons 
N.Y.  - - 
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We  now  represent  the  path  dependency  for  the  emitted  acoustics  In  a  homo¬ 
geneous  atmosphere  as 


G(r,t.  0.tQ) 


lR./c„  -  -  to)] 

R 


with  R,  (t  -  tQ)  >  0. 

Near  the  ground  surface,  V^(t)  «  2  for  all  but  near  the  air-coupled  frequen¬ 
cies,  To  these  approximations,  the  measured  surface  pressure  at  an  offset  A 
Is  taken  to  be: 

p<A,t)  *  2.H(t)*6(R/ca  -  (t  -  to))/R*P^tQ)*e<t0)  .  (4) 

For  a  vertical  launch,  the  surface  pressure  Is 

PA(t  +  A  /c2)  «  2/A(c/c  )  •  P,(t)*e(t)  (5) 

A  c1  a  '  or  <p 

o 

where  the  delay  time,  Ac^/c,  and  gain  factor,  2/A(c/0,  are  time  dependent  be¬ 
cause  c  =  c(t). 

In  turn,  pressure  amplitude  ratios  across  our  array  ore 

IV'V^AjI  =1-<ca/°)2‘6A1/^  +  6Aj) 

with  the  delay  times 

6  /c  .  (7) 

1 

Much  as  before  we  test  our  representation  of  surface  pressure  through  an 
ensemble  of  "best"  amplitude  and  delay  pairs  by  fitting  our  observations  in  a  least 
squares  sense  to  the  form  of'  Eqs.  (0)  and  (7).  The  best  amplitude  ratios  (Fig¬ 
ure  22),  and  delay  times  (F’lgure  23),  are  obtained  for  a  2.  56-sec  data  time  gate 
during  peak  surface  loading.  The  amplitude  ratios  and  delay  times  so  obtained  are 
consistent  with  an  "equivalent"  plume  source  at  an  elevation  of  350  ft.  (107  m), 
well  below  the  Titan  rocket. 

We  again  evaluate  our  representation  by  computing  the  relative  error  between 
the  predicted  and  observed  pressure  values  with  separation,  as  seen  in  Figure  24. 
The  error  row  is  a  sensitive  function  of  separation.  Indeed,  for  a  separation  of 
250  ft.  (76.2  m),  the  mean  square  error  is  about  50  percent  of  the  observed  value. 
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Figure  22.  Amplitude  Ratio  During  Peak  Loading 
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Figure  23.  Phase  Velocity  During  Peak  Loading 


Figure  24.  Residual  Error  for  Maximum  Loading 


By  Ignoring  the  motion  of  the  rocket  over  the  sample  time  gate,  we  misrepresent 
the  actual  phase  change  In  the  acoustic  load.  The  best  average  phase  shift  under¬ 
estimates  the  true  shift  at  the  beginning  of  the  sample.  Towards  the  end  of  the 
sample,  the  best  average  value  Is  too  targe. 

4.4.3  PHASE  VELOCITY 

In  order  to  sqrt  out  plume  and  propagation  properties,  we  make  a  succession 
of  phase  velocity  estimates  for  the  plume  generated  pressure  propagating  over  our 
array.  Figure  25  Is  a  set  of  median  phase  velocities  calculated  from  the  phase 
shifts  of  a  sequence  of  overlapping  2.  56-sec  samples  (Appendix  E). 

We  now  set  out  to  show  that  the  propagation  Is  nondisperslve;  that  Is,  that  the 
phase  velocity  Is  Independent  of  frequency.  As  computed  from  the  cross  spectra, 
"phase  velocity"  Is  a  set  of  magnitudes  ordered  by  frequency  (see  Figure  26).  We 
now  treat  the  calculated  velocities  to  be  Inconsistent  independent  estimates  about  a 
constant  velocity.  Under  this  assumption,  we  can  expect  that  the  mean  number  of 
runs  computed  from  sign  changes  In  the  first  differences  In  our  phase  velocity 
estimates  is  l(2n  -  1)/3|  where  n  is  the  total  number  of  estimates.  23  For  large  n 
the  number  of  runs  tends  to  be  normally  distributed  about  the  median  with  a  stand¬ 
ard  distribution  of  [(16n  -  29)/901 1/2.  For  our  tests  n  =  127,  hence  the  expected 
number  of  runs  is  84  with  a  standard  deviation  of  5.  For  nondlspersive  propagation 

23.  Owen,  D.  (1962)  Handbook  of  Statistical  Tables.  Addison-Weslej  Publlshine 
Co.,  Reading,  MA. 
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Figure  25.  Estimated  Phase  Velocity  vs.  Time 


we  can  expect  that  for  about  VO  percent  of  our  trials  the  number  of  runs  will  lie  In 
the  range  79  sfii  s  59.  Figure  27  gives  the  actual  number  of  runs  obtained  for  the 
launch  aa  they  relate  to  sample  start  time.  We  conclude  that  acoustic  phase  ve¬ 
locity  over  the  array  is  Independent  of  frequency.  Plume  generated  propagation 
delays  can  he  calculated  simply  from  <H/e  )  =  t.  Also,  the  "equivalent"  plume 
Source  height  Is  Independent  of  frequency.  Our  calculations  do  not  support  the 
notion  that  the  "source  pressure"  Is  refracted  through  an  angle  proportional  to 
wavelength. 

4.4.4  COHERENCE 

An  estimate  of  coherency  based  on  observations  taken  over  the  entire  launch 
sequence  is  calculated  after  the  manner  of  Hinlch  and  Clay*5  for  a  station  pair 
located  200  ft.  (61  m)  apart  (Figure  28).  This  ordinary  estimate  of  coherency 
leads  to  an  "Incoherent  pressure  field"  for  all  but  the  lowest  frequencies.  The 
low  coherency  values  found  here  in  our  primary  measurement  band  (3-30  Hz) 
reflect  the  nonstationarlty  of  the  amplitude  and  phase  relation  of  the  surface  pres¬ 
sure  during  launch.  To  show  that  this  is  indeed  the  case,  we  now  calculate  the 


"spacial  coherency"  of  the  surface  pressure  for  the  same  station  pair  by  compen¬ 
sating  our  pressure  phase  measurements  In  accordance  with  our  estimate  of  c(t) 
and  the  relations  of  Eq.  (5)  and  (7).  Figure  29  is  our  estimate  of  coherency  after 
the  pressure  measurements  are  weighted  and  shifted  to  agree  with  our  estimate 
of  c(t).  The  pressure  field  la  now  quite  coherent  up  to  30-35  Hz.  A  stationary 
amplitude  and  phase  relation  results  between  field  points  in  our  primary  passband 
when  the  data  are  corrected  for  rocket  climb-out. 

The  loss  of  coherency  above  35  Hz  in  Figure  29  is  readily  accounted  for  by 
the  ambient  r.oise  measurement  and  the  low  protection  ratio  for  the  digitized  data 
much  above  3  5  Hz.  As  can  be  seen  In  Figure  30,  the  upper  limb  of  the  plume  gen¬ 
erated  acoustic  load  falls  into  the  ambient  noise  shortly  after  the  time  of  peak 
loading. 

In  our  passband  we  find  that  the  plume  generated  acoustic  load  is  spaclally 
coherent  for  stations  separated  by  as  much  as  200  ft.  (til  m). 

4.4.5  SPECTRAL  ESTIMATES 

Stable,  high  degree  of  freedom  (DOF)  plume  spectral  estimates  at  some  ele¬ 
vation  angle  $  can  be  obtained  in  a  number  of  ways.  One  obvious  technique  is  to 
Increase  the  DOF  and  thus  acquire  stability  by  averaging  spectra  over  a  number 
of  launches;  that  is,  ensemble  average.  In  turn,  for  any  one  launch  an  estimate 
can  be  stabilized  by  averaging  over  an  ensemble  of  measurements.  Figure  31  is 
our  spectral  estimate  during  peak  load  obtained  by  averaging  perlodograms  over 
eight  stations. 

Lastly,  we  can  increase  DOF  by  averaging  over  frequency.  Here  we  attempt 

to  construct  high  DOF  spectra  by  fitting  perlodograms  to  a  spectral  form  over  our 

24 

entire  measurement  passband.  Both  from  physical  considerations  and  a  number 

1 9 

of  experimental  studies,  the  power  spectral  density  of  the  surface  pressure  dis¬ 
turbance  emitted  by  an  undeflected  chemical  rocket  plume  is  found  to  closely 
follow  the  form 

P  (id)  ■  4  P  /ffld  (ij  /  Id  +  Id  /id)  ^  (8) 

pp  o'  m  '  m  nr 

with  P0/7rwnl  the  value  at  the  spectral  maximum,  um  the  angular  frei  ency  at  the 
spectral  maximum,  and  Pq^  the  overall  sound  power. 

We  now  test  the  hypothesis  that  the  acoustic  power  spec  tra  density  emitted  t  ■ 
the  steady-state  undeflected  Titan  ill-D  plume  can  be  east  in  the  form  equivalent 
to  Eq.  (8).  The  pressure  measured  at  a  surface  station  950  ft.  (290  in)  from  the 


24.  Powell,  A.  (1054)  Theory  of  vortex  sourd,  J.  Aeoust.  Soc.  Am.  35:177-195. 
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launch  pad  is  used  to  generate  a  succession  of  two  DOF  periodogram  estimates 

based  on  measurements  taken  in  2.  56-sec  sample  gates.  The  spectral  estimates 

calculated  here  are  unbiased,  inconsistent  estimates  of  a  stationary  stochastic 

process  when  the  time  samples  are  concatenated  with  a  like  number  of  zeros 
2  5 

before  transforming.  The  spectra  constructed  here  are  power  spectra  in  the 
sense  that  they  satisfy  Parseval's  relation  in  the  time  interval  of  the  sample, 
namely 

<e2>T  ■/  V')d,*mo 

o 

where  T  =  2.  56  sec  and  f  is  the  Nyquist  frequency. 

The  spectral  form  of  Eq,  (8)  <s  fitted  to  each  periodogram  to  obtain  a  best 

estimate  of  (Pq^  fm)  in  the  sense  that  the  mean  square  difference  between  the 

two  representations  is  minimized  by  our  choice  of  Pn . ,  f  .  Figure  32  shows  the 

uii  m 

best  spectra  calculated  around  the  time  of  maximum  acoustic  loading,  superim¬ 
posed  on  the  calculated  periodogram  points. 

We  now  construct  a  table  of  (P^,  f  )  values  over  a  sequence  of  overlapping 
2.  53  sec  gated  samples  (see  Table  1).  Also  included  in  Table  1  is  the  overall 
mean  square  surface  pressure  estimate  developed  by  the  model  rocket  plume.  It 
is  worth  not.ng  here  that  the  theoretical  total  pressure  level  is  less  than  the  meas¬ 
ured  value  for  sample  midtimes  after  15  sec.  At  such  times  the  acoustics  are 
abating  while  the  turbulence  at  the  site  is  well  developed. 

A  number  of  tests  were  conducted  to  evaluate  the  hypothesis  that  the  meas¬ 
ured  spectra  did  not  depart  significantly  from  the  form  of  Eq.  (8).  One  such  test 
result  is  given  In  Figure  33.  The  figure  is  a  probability  paper  plot  of  the  distribu- 
tion  of  the  observed  residuals  against  the  residual  distribution  expected  for  a 
variate  with  two  DOF:  the  expected  distribution,  had  we  fitted  the  periodogram 
ordinates  to  the  true  spectra.  The  paper  is  constructed  so  that  the  acceptance 
criteria  is  one  of  simply  accepting  the  plotted  values  to  lie  on  the  Indicated 
straight  line.  Both  from  such  tests  and  tests  based  on  an  acceptance  range  in  our 
Figure  of  Merit  obtained  from  simulation,  we  conclude  that  the  periodograms 
are  well  represented  by  an  equation  of  the  form  of  Eq.  (8)  particularly  around  the 
time  of  maximum  acoustic  loading. 


25.  Ywen,  C.  (1979)  On  the  smooth  periodogram  method  for  spectrum  estimation, 

Signal  Processing  1:83-86. 

26.  Hartnett,  E. ,  and  Carleen,  E.  (1980)  Characterization  of  Titan  III-D  Acoustic 

Pressure  Spectra  by  Least  Squares  Fit  to  Theoretical  Model.  AFGL-TR- 

‘80-004,  AD - -  - 
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Table  1.  Overall  Sound  Power  Values 


Start  Time 
After  Ignition 
(see) 

f 

o 

(list) 

— 

Th«*orntlcal 
P0(pH'2>  db 

Error 

M) 

Observed 

Power 

(psl)2 

o  <  f  <  50  Hz 

Spectral 

Maximum 

Figure 

of 

Merit 

3.  83 

(i.  708 

8.477(10”°)  130,0 

4.022U0'10) 

7.931(10'°) 

214.  8 

5.30 

14.200 

(i.  741110'4)  130.  0 

1.2GB(10'a) 

3.875(10'°) 

1.  5008(10'°) 

1.005 

(>«  95 

12. 071 

1.  00U(10':1)  142.  8 

5.342<10*8) 

1.007U0'3) 

4.  2350(10'°) 

0.  8204 

8,  51 

11. 100 

1.  257(10'3)  141.7 

0.848U0'8) 

8.284U0'4) 

3.  578(10'°) 

1.381 

10.  07 

7.  548 

8.  30l(10'4)  138.  8 

1.  U48(10”8) 

4.  003(10'4) 

2.  883(10'°) 

1.  528 

11.  (13 

0.  235 

2.  845(10"4)  135.3 

S.418(10'9) 

2.  572(104) 

1.730(10'°) 

1.  813 

13.  10 

(1.  007 

1.  116(10'4)  131.  2 

7.  147(10”10) 

1.  035(10'4) 

5.  821(10'°) 

1.  405 

14.  75 

4.  300 

4.  UOB<10'°)  127,4 

2.  108(10~10) 

4. 872(10'°) 

3.341(10'°) 

2.  50 

1(1,31 

. 

2.855 

1.  8(I0<  10-B)  123.4 

- 

0.  203(10'U) 

2.420(10'°) 

2.  074(10'“) 

1.  032 

We  now  plot  PQA  as  it  relates  to  the  start  time  of  the  sample  following  rocket 
Ignition  In  Figure  34.  The  time  for  a  maximum  PqA  f or  an  observer  located 
£>50  ft  (200  m)  from  the  launch  pad  Is  about  8  sec.  Following  this  maximum,  the 
overall  sound  power  estimate  decays  exponentially  at  a  rate  of  about  2.  5  db/sec. 

In  Figure  35,  fm  is  plotted  as  it  relates  to  the  sample  start  time.  Over  the 
period  of  estimation,  the  frequency  of  the  spectral  muxlmum  decreases  roughly 
linearly  with  time.  The  one  sigma  spread  In  fm,  found  by  simulation,  Is  Indicated. 
The  shift  In  t'm  Is  significant.  The  plume  emits  significantly  lower  frequency 
spectra  at  the  higher  elevation  angles.  The  shift  is  well  in  excess  of  the  frequency 
shift  ascribable  to  the  source  velocity. 


_ 1 _ I _ L, 

9  10  19 

SAMPLE  START  TIME  (MCOndt) 

Figure  34.  Overall  Sound  Power  Level 


TIME  AFTER  IGNITION 

Figure  35.  Frequency  at  the  Spectral  Maximum 


5.  ADMITTANCE  ESTIMATES 


Over  the  first  20  sec  of  launch,  admittance  estimates,  based  on  averaging 
three  consecutive  0.  64-sec  co-located  pressure  and  ground  samples,  are  com¬ 
puted.  The  computation  is  repeated  every  1.  56  sec  following  the  arrival  of  the 
ignition  pulse.  As  in  the  theoretical  case,  we  again  locate  admittance  maxima  in 
terms  of  phase  velocity  and  frequency.  We  find  that  the  response  maxima  mi¬ 
grates  in  the  (c,  f)  plane  in  the  expected  manner  in  that  the  "frequency"  of  the  air- 
coupled  disturbance  decreases  with  an  increase  in  phase  velocity  for  both  the 
fundamental  and  the  first  higher  mode  response,  shown  in  Figure  3  6.  In  turn,  the 
absolute  value  of  the  ground  response  maxima  are  given  as  they  relate  to  phase 
velocity  in  Figure  37.  The  measured  responses  are  substantially  less  at  the 
higher  phase  velocities,  as  expected. 

We  now  calculate  the  usual  plume  directivity  value*3  based  on  our  measured 
overall  plume  pressure  levels,  PQ^,  and  1/R  attenuation  for  a  source  located 
near  the  base  of  the  Titan  (see  Figure  36).  If  we  now  compare  this  pattern  to  the 
radiation  pattern  obtained  using  Eq,  (5),  an  elevation  angle  obtained  from  our 
estimate  of  phase  velocity,  we  find  the  usual  directivity  calculation  poorly  repre¬ 
sents  the  pressure  load  for  facility  response  computations  (Figure  36). 


Figure  36a.  Measured  Admittance:  Fundamental 
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Figure  3Gb.  Measured  Admittance:  First  Mode 


Figure  37b.  Admittance  Maxima:  First  Mode 
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Figure  38. 
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Figure  39.  Radiation  Pattern  From  Phase  Velocity 

6.  SCALING  TITAN  US  D  INTO  STS 
6.1  Power  S[xctni  Demit y  E*tim»te» 

The  overall  acoustic  efficiency  of  undo  fleeted  chemical  rocket  plumes  is  rela¬ 
tively  constant  over  a  wide  range  of  thrust  levels  for  similar  nozzle  flow  param¬ 
eters.  Table  2  summarizes  the  engine  parameters  needed  to  scale  the  acoustic 

21  27 

loads  due  to  Titan  II1-D  Into  STS  launches.  '  We  now  normalize  our  estimate 
of  the  Tital  II1-D  spectra  observed  at  950  ft.  (290  m)  around  peak  loading  using  the 
scaling  parameters  given  In  Table  2.  In  standard  form,  shown  In  Figure  40,  the 
maximum  spectral  value  is  expected  at  f(de/X.Je)  =  0.02;  our  measured  value  is 
0,0198.  In  standard  form  we  expect  the  spectral  maximum  W(f)/WOA  •  (Ue/de)  to 
be  10  db  in  excess  of  the  overall  power.  Our  measured  value  is  12  db.  Our 

27.  Smith,  J. ,  The  Aerospace  Corporation,  Ltr  dated  4  January  1980,  Rocket 
Engine  Parameters. 


Table  2.  Rocket  Engine  Parameters 


Number 

of 

Engines 

Nozzle  Exit 
Diameter  (de) 
in.  (cm) 

Exhaust  Velocity 
at  Sea  Level  (Ue) 
ft/sec  (m/eec) 

Total  Power 
at  Liftoff 
lb  ft/sec 
(kg  m/sec) 

STROUHAL# 

SRB 

2 

141.7(359.  9) 

6440(1963) 

3.  3 1710(4,  5969) 

- 

SSME 

3 

92.0(234) 

0730(2966) 

1.  21810(1.  6879) 

- 

STS 

- 

- 

- 

4.  53510(6.  2839) 

0.  00293f 

TIH-D 

2 

106.  6(270.  8) 

7676(2340) 

1.  92110(2.  6619) 

0,  001637f 

spectra  for  a  Titan III-D  satisfies  the  standard  form  spectra  for  undeflected 
chemical  rocket  plumes  in  the  thrust  range  350  to  7,  000,  000  lb  (159  to 
3,  180,  000  kg).20 

We  now  construct  a  normal  stationary  time  series  that  fits  the  spectral 
properties  of  the  surface  pressure  expected  at  a  point  950  ft  (290  m)  from  the  pad 
for  a  STS  launch  around  the  time  of  maximum  loading.  The  design  OASPL  for  the 
STS  Simulation  is  2.  36  times  that  of  the  Titan  III-D;  that  is,  2.36  is  the  ratio  or 
the  total  power  levels  between  the  two  systems  at  launch.  In  Figure  41  we  com¬ 
pare  the  spectra  for  the  two  systemB.  Each  spectra  satisfies  the  standard  form 
spectra  shown  in  Figure  40,  when  the  scaling  values  of  Table  2  are  applied.  We 
find  that  for  structural  responses  above  20  Hz  the  expected  motion  levels  due  to 
STS  and  Titan  launches  differ  only  by  a  few  db.  For  low  mode  ground  and  large 
structure  responses,  the  expected  motion  environments  of  the  two  Bystems  differ 
significantly;  that  is,  <20  Hz. 

In  order  to  address  the  problem  of  repeated  launches  we  construct  STS  spec¬ 
trum  estimates  at  the  50  and  97.  5  percent  confidence  levels  by  appealing  to  the 

2  9  G 

ergodlc  hypothesis  and  the  x  distribution  of  our  test  statistics. 

X2(x)  =  F(2Pj/P(fj))  DOF  *  2 

For  example,  the  lower  bound  spectra  PB0(f|)  shown  in  Figure  42  will  be 
exceeded  50  percent  of  the  time,  on  the  average.'  by  a  perlodogram  estimate  at 
the  frequency,  f,  when  ensemble  averages  equal  time  averages, 

6.2  OBSPL  Eitimatci 

We  now  compare  our  forecast,  bn3ed  on  scaling  up  power  spectra  obtulned 

from  one  Titan  launch  around  the  time  of  maximum  acoustic  loading,  with  OBSPL 

estimates)  obtained  by  scaling  down  ensemble  averages  of  Saturn  V  launches  during 
2 1 

the  liftoff  peak,  In  order  to  compare  the  two  "spectra,  "  a  number  of  adjustment 
must  be  introduced.  The  factors  needed  to  convert  our  results  Into  OBSPL  for 
different  confidence  levels  are  given  in  Table  3. 

The  table  y^  quantities  given  here  are  needed  to  adjust  the  x2  distributions 
of  the  OBSPL  estimates  when  Integrating  constant  resolution  spectra;  that  Is, 


DOF-Yj-  *  -dof72 - *  f  x(DCil''2)/?  •  e<_X/2>  dx  «N*/100 

m  2DOF/z  .  s/DOF/2  J 


where  50,  97.  5. 
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Figure  4,'i.  STS  Surface  Spectra 


Table  3.  Overall  3ound  Pressure  Levels 


fc 

E{OBSPL}(2.  56) 

DOF 

V97. 5% 

2 

129.4 

6 

0.892 

2.408 

4 

137.2 

12 

0.942 

1.945 

8 

141.3 

26 

0.973 

1.612 

16 

141.7 

54 

0.987 

1.411 

315 

140.2 

106 

0.993 

1.287 

63.0 

137.  5 

214 

0.997 

1.  189 

Figure  43  summarizes  our  STS  OBSPL  forecast  for  surface  pressure  at  the 
50  percent  and  97.  5  percent  confidence  levels.  In  the  usual  fashion,  the  dis¬ 
persion  about  the  expected  value  Is  sensitive  to  the  DOF  of  the  estimate.  For  an 
averaging  time  as  large  as  2.  56  sec,  the  expected  and  median  spectral  estimates 
differ  by  less  than  5  percent  above  8  Hz.  For  bandwldths  somewhat  larger  than 
30  Hz,  there  Is,  on  the  average,  only  a  1  db  difference  batween  the  50  and  97.  5 
percentiles  for  the  two  distributions. 

We  now  must  allow  for  differences  In  our  OBSPL  due  to  differences  In  averag¬ 
ing  times,  In  Figure  44  we  plot  the  maximum  of  the  RMS  pressure  level  found 
over  passband  for  the  Titan  iaunch  as  It  relates  to  averaging  time.  For  the  Titan, 
we  have  a  3.3  db  separation  between  the  RMS  pressure  computed  for  a  0.3  sec 
averaging  time  and  one  based  on  2.  56  sec.  For  large  DOF  estimates  (greater 
than  25)  we  now  simply  add  3.3  db  to  our  results  given  In  Figure  43  to  compensate 
for  time  averaging  "on  the  order  of  0.3  sec"  over  a  lift  of  peak  that  lasts  for 
2.3  sec  (based  on  3  db  points  of  the  SATURN  V  OASPL  envelope),21  Our  results 
now  are  quite  close  to  the  NASA  mean  OBSPL  curve  for  STS  obtained  by  down 
scaling  SATURN  V  data.  All  our  results,  seen  In  Figure  45,  are  well  below  the 
NASA  97.7  percent  confidence  level  forecast.  Indeed,  the  constant  3.  1  db  separa¬ 
tion  between  the  50  and  97.7  percent  confidence  level  appearing  iere  cannot  be 
explained  In  terms  of  the  equivalence  of  ensemble  and  time  avernges. 


48 


160 


Figure  45.  STS  Estimates  OBSPLi  <  )  0.  3  sec 


7.  GROUND  RESPONSE  AT  SLC  6 

The  ground  response  in  the  vicinity  of  the  Launch  Control  Center  (LCC)  at 
SLC-6  wa3  measured  for  an  acoustic  source  located  300  ft  (91.4  m)  above  the 
STS  launch  pad,  as  shown  in  Figure  46.  The  ground  admittance  measured  over 
this  path  is  dominated  by  the  fundamental  air-coupled  wave.  The  frequency  of  the 
air-coupled  term  found  near  the  LCC  is  at  6  Hz.  This  frequency  is  significantly 
higher  than  the  frequency  for  the  fundamental  air-coupled  wave  excited  at  SLC-4. 
The  ready  explanation  for  the  difference  lies  in  a  difference  in  the  thickness  of 
the  unconsolidated  sediments  at  the  two  sites. 

We  now  construct  a  sample  STS  pressure  P(t)  to  simulate  the  expected  maxi¬ 
mum  loading  for  an  elevated  plume  source  for  the  distance  to  the  LCC;  that  is, 


p(ALCC,t)  =  2/4LCC(c/ca)  •  P^(t>  >1'  £(t) 


where  as  before  2(t)  is  from  a  N(o,  a),  stationary  process  and  (  P^  *  P^)  satisfies 
the  standard  form  spectral  figure  given  in  Figure  41  for  the  STS  at  altitude. 

The  sample  pressure  in  turn  is  passed  through  an  operator  that  satisfies  the 
measured,  path  sensitive  admittance  to  give 


LCC 


(t)  =  Y, 


(t)  *  P, 
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LCC 
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location  a 


Figure  4U.  Ground  Admittance  near  LCC,  SLC-6 


the  predicted  ground  motion.  The  resulting  "ground  motion  samples"  can  then  be 
used  to  characterize  the  STS  motion  environment  around  the  time  of  maximum 
loading. 

The  elastic  vibrations  expected  within  the  LCC  due  to  launch  generated 
acoustics  can  be  forecast  in  much  the  same  manner  that  we  have  just  constructed 
the  ground  term.  As  in  the  case  of  ground  motion,  the  key  element  is  the  correct 
representation  of  the  spacial  and  temporal  properties  of  the  acoustic  load.  A 
test  explosion  located  above  the  STS  launch  pad  ..as  used  to  calculate  the  "impulse 
response"  of  a  point  within  the  structure  when  the  test  and  launch  generated 
acoustics  have  common  loading  characteristics  over  the  structure  (see  Figure  47). 
The  vibration  excited  during  a  STS  launch  is  then  estimated  by  temporal  convolu¬ 
tion  with  the  acoustic  load  anticipated  for  a  STS  launch. 

The  effect  of  the  LCC  on  ground  level  motion  can  be  seen  by  forming  the 
ratio  between  the  LCC  motion  amplitudes  found  at  ground  level  with  correspond¬ 
ing  amplitudes  obtained  for  the  ground  outside  near  the  LCC  (see  Figure  48). 

From  Figure  46  it  is  clear  that  the  exterior  ground  motion  is  dominated  by  an 
air-coupled  wave  near  5-8  Ilz.  In  this  band  the  corresponding  amplitudes  within 
the  LCC  at  ground  level  are  less  than  the  exterior  terms.  The  structure  has  the 


effect  of  attenuating  the  air-coupled  wave.  In  contrast  over  the  band  12-28  Hz, 
the  building  responses  clearly  dominate.  Motion  in  this  band  is  substantially 
larger  than  the  corresponding  ground  levels  computed  outside. 

Where  comparisons  are  possible,  motion  estimates  based  on  linear  elastic, 
three-dimensional  finite  element  modeling  of  the  LCC  moderately  overestimates 
both  the  level  and  duration  of  the  basement  motion  due  to  seismo-acoustic 
loads.  The  damping  terms  selected  for  modeling  the  LCC  response,  while 
within  the  range  of  accepted  values,  are  found  to  be  conservative.  It  is  quite 
pruder'  that  damping  terms  be  selected  conservatively,  for  structural  response 
to  seismic  excitation  is  quite  properly  a  major  subject  area  for  research. 
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29.  Simpson,  1.  (1978)  On  the  interaction  of  Rayleigh  waves  with  structures, 
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List  of  Symbols  and  Abbreviations 


c 


e  ,  U 
o' 

DOF 

f 

fm 

fn 

FFT 

FFT- 

G 

H 

Hs 

Hn 

H(t) 

M 


m 

o 

n 

OASPL 

OBSPL 


phase  velocity 

phase  velocity:  air  wave 

phase  velocity:  ltuylelgh  wave  at  array  offset  distance 
phase-group  velocity  for  model  half  space 
degrees  of  freedom 
frequency 

maximum  spectral  frequency 

Nyqulst  frequency 

Fast  Fourier  Transform  (forward) 

Fast  Fourier  Transform  (Inverse) 

Green's  function 
depth 

source  height 
rocket  height 

measurement  system  response 

figure  of  merit 

mean 

number  of  runs 

overall  sound  pressure  level 

octave  band  sound  pressure  level 

overall  sound  power 

pressure  power  spectral  density 

source  shaping  operation 
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pa 

p' 

Pi 


Pr 


P 

Vs 

Vz 

Vz 

V.(t) 


w(r) 

WO  A 
Y(c,  f) 

Z(c,  f,  z  =  o) 
>l< 

(  )  Sec 
A 
Aj 


6  A. 


'j 


j  v 

•<‘o> 


pressure  amplitude 
predicted  pressure 
Incident  pressure 
reflected  pressure 
surface  pressure 
source  observer  distance 
source  distance  from  origin 

mean  square-reference  pressure 

duration 

time 

reference  time 

ith  sample  time 

compressloned  wave  velocity 

shear  wave  velocity 

particle  velocity 

particle  velocity  estimate 

ground  surface  pressure  operator 

exhaust  power 

overall  exhaust  power 

acceptam 

impedance 

linear  convolution 

time  averaging 

offset  source  distance 

offset  launch  pad  distance 

reference  launch  pad  distance 

Dirac  delta  function 

reference  point  offset 

mean  square  error  -  jth 

zero  mean,  normal,  stationary,  independent  stochastic  process 
density 

rms  of  normal  process 
delay  time 
best  delay 

observer-source  elevation  angle 
angular  frequency 


Appendix  A 

Computerized  Algorithm  for  a  Layered  Whole  Space 

The  method  of  ■Abo~Zena,)  was  used  to  compute  the  behavior  of  acoustic  waves 
In  a  fluid  half-space  (air)  Incident  on  an  alluvial  layer  resting  on  a  half-space  of 
rock  (Figure  1). 

For  each  of  the  media: 

y  =  -Ikb^  cosh  (r^kz)  -  Ikb^  slnh  (r^kz)  -  7^kbg  slnh  (r^kz) 

-  r^kb^  cosh  (r^kz) 

w  =  r  kb,  slnh  (r  kz)  +  r  kb„  cosh  (r  kz)  -  lkb„  co3h  (r„kz) 
o'  1  a  a  &  a  o  p 

-  Ikb^  slnh  (r^kz) 

a7  =  p(y  -  1)  k2c2b^  cosh  (r  kz)  +  f)(y  -  1)  k2c2  sinli  (r^kz) 

e>  2 

-  Ik-'r^pbg  slnh  (r^kz)  -  ir^2ph  cosh  (r^kz) 

T  =  -2lp(32r^k2b1  slnh  (r^kz)  -  2ir j)ft^k2h2  cosh  (r^kz)  -  (7  -  l)pk2c2bg  cos 


where 


u  «x  component  of  displacement 

w  3  z  compom.it  of  displacement 

a  «  normal,  stress 

z 

T  “  tan8enfclal  stress 

p  -  density 

=  Lam <5  elastic  constants 


c  a  phase  velocity  of  free  wave  along  x  axis 
w  =  radial  frequency 


>■,3 

-  ’Ml)2 


and  bp  b2,  and  are  to  be  determined  for  each  medium  according  to  boundary 
cond  it  Iona. 

Using  Abo- Zena's  method  and  notation,  the  system  of  equations  for  each 
layer  can  be  combined! 

B„  =  U.,(0)-1A.,Ur»<0)B„ 


^ '  Pt«’!\rr*-7  *  ?'  'T/?  nr<  '-TO\vrw  vr>i  «■;■  ;w\'l  rr  «-  '■M-(^M,.;.,„l^„„r.,ir. 


■■.{•■’ r'  l  M'.«-;  t  v  ■  * 


?.*nwnT?(^|i'"‘>-'  »>  v ■ 


where 


=>  set  of  b's  for  top  layer  (air) 


■  set  of  b's  for  bottom  layer  (rock) 


V1'  0  0  (,D^0r30 

0  P0y0r°o  p0%_1)  0 


U2(0)  r. 


-(72-D 


(7g  -  1) 


‘>2% 


A,  ■  vAhWAO  V 


where 


00 


and  since  this  layer  Is  semi -Infinite  bg  =  b^  =  0  multiplying: 


G2 


where 

'  V2 

*1  ■  -P0b! 


and  x  is  an  indefinite  term  due  to  the  »  •  0  product.  Representing  the  U„(0)_1A 

z 

product  as 


11 

c12 

c13 

c14 

21 

c22 

C23 

c24 

31 

c32 

c33 

c34 

41 

c42 

c43 

c44 

the  equations  can  be  rewritten! 

B2  *  CY  * 

In  the  bottom  medium,  to  ensure  waves  traveling  in  the  positive  z  direction  (no 
source  at  z  -  »),  b^  =  and  b|  =  b*  therefore, 


bl 

X 

bl 

. 

C 

, 

y2 

-jbg 

yl 

-ibl3 

0 

-  — 

M,  — 

_  _ 

x  can  be  eliminated  and  Solve^  for: 


y2/yi 


(n41  "  c31'|c13  ~  c23|  +  |c43  "  c33)|c21  "  cll| 
^c31  ”  °41^c12  '  c22)  +  (c32  “  c42^c21  "  cll* 


and 


b2/bl  =  ”P0 /r«Q  •  Vyl 
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Appendix  B 

Inverts  Operator 


System  response  was  fine  tuned  by  generating  theoretical  responses,  com¬ 
puting  inverse  FFT's  and  comparing  these  with  calibration  responses.  Overall 
gain,  pole  location  of  seismometers,  and  damping  were  adjusted  until  the  differ¬ 
ence  between  theoretical  and  observed  was  minimized  in  the  least  squares  sense. 

Having  generated  the  theoretical  response  of  the  system,  the  next  task  was 
to  generate  an  operator  which  would  yield  the  input  of  the  system  when  convolved 
with  its  output.  To  do  this,  the  reciprocal  of  the  magnitude  of  the  response  in 
the  frequency  domain  was  computed  along  with  the  negative  of  the  phase.  To 
Insure  that  the  Inverse  FFT  would  be  real,  a  linear  function  of  frequency  was 
added  to  the  phase  such  that  the  phase  at  the  Nyqulst  frequency  would  be  -180  de¬ 
grees.  To  reduce  the  effect  of  corners,  the  value  of  the  magnitude  at  the  Nyqulst 
frequency  is  halved.  The  inverse  FFT  was  computed  using  512  samples.  The 
operator  thus  produced  had  enough  zeroes  to  ensure  that  the  convolution  could  be 
performed  linearly. 

To  check  the  performance  of  the  operator,  a  forward  operator  was  generated 
by  computing  the  inverse  FFT  of  the  system  response.  This  was  convolved  with  a 
series  of  random  numbers.  When  the  result  was  convolved  with  the  inverse 
operator,  the  series  produced  had  an  RMS  within  0.  1  percent  of  the  original. 


Appendix  C 

Admittance  Estimate  and  Prediction  Operator 


To  estimate  the  admittance,  the  pressure  and  vertical  ground  velocity  at  the 
arrival  of  the  Ignition  pulse  were  Isolated.  Using  speed  of  sound  in  air,  the 
estimated  arrival  time  was  0.83  sec  after  ignition.  The  corresponding  channels, 
corrected  for  system  response,  were  multiplied  by  the  following  factors 

“0 

1  -  cos  (t  -  3)/4  •  ff 

- '-2  - 

f  ■=  1 

1  +  cos  (t  -  2.  l)/0.  4  •  it 

-  2~.- . . 

0 

FFT's  containing  512  points  were  comouted,  and  the  ratio  of  amplitude  of  velocity 
to  amplitude  of  pressure  and  difference  between  phase  of  velocity  and  phase  of 
pressure  used  to  generate  an  operator  by  computing  the  inverse  FFT.  This 
operator  continued  nonzero  throughout  its  5.  12  sec  length.  To  eliminate  this, 
it  was  multiplied  by  a  decaying  exponential  thus: 


t  <  0.  3 
0.  3  <  t  <  0.  7 
0.  7  s  t  <  2.  1 
2.  1  <  t  <  2.  5 
t  >  2.  5 
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Appendix  D 

Titan  III  Launch  Profile 
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Table  Dl.  Titan  111  Launch  Profile  Values 


Simulated 

t(sec) 

Range  In  ft.  (m) 

log^t 

log10R 

1.015 

5(1.5) 

0.0065 

0. 6990 

1.404 

12(3.7) 

0.1474 

1.0792 

1.795 

22(6.  7) 

0.25^1 

1.3424 

2. 185 

34(10.4) 

0.3395 

1.5315 

2.  575 

50(15.2) 

0.4108 

1.6990 

2.966 

68(20.  7) 

0.4722 

1.8325 

3.356 

89(27.  1) 

0.5258 

1.9494 

3.747 

113(34.4) 

0.5737 

2.0531 

4. 137 

140(42.7) 

0.6167 

2.1461 

4.527 

170(51.8) 

0. 6558 

2.2304 

4.918 

202(61.  6) 

0. 6918 

2.3054 

5.308 

238(72.  5) 

0.7249 

2.3766 

6.089 

319(97.2) 

0.7845 

2.5038 

6,870 

413(125.9) 

0.8370 

2.6160 

7.  651 

520(158.  5) 

0.8837 

2.7160 

8.432 

639(194.  8) 

0.9259 

2.8055 

9.212 

771(235) 

0.9644 

2.8871 

9.  993 

916(279.2) 

0.9997 

2.9619 

10.774 

1075(327.7) 

1.0324 

3.0314 

12.336 

1433(436.  8) 

1.0912 

3.1562 

13. 116 

13.897 

14. 678 

15.459 

2316(705.9) 

1.1892 

3.3647 

16.240 

17.021 

17.801 

3127(953.  1) 

1. 2o04 

3.4951 

18.582 

18.973 

20. 144 

4071(1241) 

1.3041 

3. 6097 

Appendix  E 

Phut  Velocity  end  Azimuth  Estimate 


To  compute  phase  velocity,  the  phase  lags  between  channels  were  computed 
by  means  of  the  FFT.  To  eliminate  ambiguities  of  2*  radians,  peak  values  of 
each  channel  were  examined  around  the  time  of  the  arrival  of  the  ignition  pulse. 
Samples  were  dropped  from  the  beginning  of  each  channel,  except  channel  15 
which  was  closest  to  the  launch  point,  until  all  the  peaks  were  aligned  and  the 
number  of  samples  dropped. from  each  recorded.  Channels  were  grouped  accord¬ 
ing  to  distance  from  each  other  along  the  x  axis  (increasing  In  west  direction). 

n 

d  (relative  number 

Channel  Pairs  x  separation  in  ft.  (m)  of  samples  dropped) 


14-11 

20(6. 1) 

2 

11-9 

20(6.1) 

2 

8-14 

40(12.  2) 

4 

9-16 

40(12.  2) 

4 

15-8 

80(24.  4) 

7 

16-13 

80(24.  4) 

7 

10-11 

0(20(12.  2)  in  y  direction) 

0 

11-12 

0(20(12,  2)  in  y  direction) 

0 

Cross  spectra  were  computed  using  FFT  and  averaged  according  to  pairing 
of  like  distances. 


X1<V  ■  l/2(F14<VFIl'V  +  Fn<VF9<V’ 

W  '  '/2<F8<VF>k>  +  WP*6<fk» 


73 


*3<V  -  »/»<p16<VF8<V  +  F16<VF»<(k» 

*y  ■  '/2iFio<yFn(v + pn«k>p*s<v» 

where  for  the  256-point  transform  of  data  sampled  at  100  samples/sec 
fk  =  (k  -  1/256)  100  . 

The  phase  angles  were  computed  and  the  phase  lag  due  to  the  original  shifts 
of  n  samples  restored: 


*j(fk)  =  tan“1(ImX1(fk)/ReX1(ffc»  +  O.Oln^*^ 

W  =  tan“1(ImX2<fk)/ReX2(fk))  +  °-01n22,rfk 
*3<fk>  =  tan_1(ImX3(fk)/ReX3(fk))  +  °*01n32,rfk 
tyfk>  “  tan_1<imY<fk)/ReY(fk)) 
where  n^  =  2,  =  4  and  ng  =  7. 

Phase  lags  were  averaged  for  a  normalized  distance  of  20  ft.  (6. 1  m)  by 
*<fk>  =  l/3(*(y  +  *(fk)/2  +  *(fk)/4) 
then: 

♦«<y  ■  [♦i«k>t*J<ik»]l/a 

and: 


V(t‘k)  =  27rfk  20/^(fk)  =  phase  velocity 
0(fk)  -  tan-1(i/iy(fk)/^x(fk))  =  azimuth  angle  . 

V  and  6  are  computed  for  each  2.  56 -sec  segment  after  ignition.  As  phase 
velocity  increases,  the  possibility  of  ambiguity  in  the  phase  lag  occurs  again. 
Phase  angles  are  monitored,  and  when  they  again  appear  to  pass  through  2ff 
radians,  the  number  of  samples  dropped  is  reduced  and  n^,  and  n^  updated 
accordingly. 


i 
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